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Abstract 



Spectra of density perturbations produced during chaotic inflation are calcu- 
lated, taking both adiabatic and isocurvature modes into account in a class 
of scalar-tensor theories of gravity in which the dilaton is metrically coupled. 
Comparing the predicted spectrum of the cosmic microwave background ra- 
diation anisotropics with the one observed by the COBE-DMR we calculate 
constraints on the parameters of these theories, which turn out to be stronger 
by an order-of-magnitude than those obtained from post-Newtonian experi- 
ments. 
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I. INTRODUCTION 



Inflationary cosmology gives a natural explanation for the horizon, flatness, homogeneity, 
and monopole problems . In its simplest form, inflation with a single scalar field predicts 
an approximately scale-invariant spectrum of Gaussian adiabatic density perturbations 0. 
If other fields are present, however, primordially isocurvature perturbations may also be 
generated to leave an imprint on the universe today 0,3- 

Scalar-tensor theories of gravitation have recently been received a renewed interest. The 
main reason is that the low-energy effective action of superstring theories generally in- 
volves a dilaton coupled to the Ricci curvature. Moreover, as pointed out by Damour and 
Esposite-Farese |^ in the study of neutron star models, a wide class of scalar-tensor theories 
not only pass the present weak-field gravitational tests but also exhibit nonperturbative 
strong-field deviations away from general relativity. Detectability of scalar gravitational 
waves by laser interferometric gravitational wave observatories(LIGO) has been explored in 



Cosmological consequences of such theories, however, have not been fully examined. 
Recently, one of us studied the generality of chaotic inflation in scalar-tensor theories in 
detail and found that the onset of inflation could be greatly affected by the presence of the 
dilaton in some cases 0. According to the attractor mechanism of scalar-tensor theories 
TI| , which works only for those theories in which the derivative of the Brans-Dicke coupling 



is positive, even those theories that are quite different from general relativity in the early 
stage of the universe approach general relativity during the inflationary stage and the matter- 
dominated stage of the universe. Even if this is the case, the deviations in the earlier stage 
may be imprinted in the spectrum of the cosmic microwave background radiation (CMB). 
We may constrain model parameters using the CMB anisotropy as a probe of the early 
universe. In fact, because of the presence of two scalar fields, i.e., the inflaton and the 
dilaton, density perturbations have both adiabatic and isocurvature modes and hence the 
spectrum may be changed. 

In this paper, we study formation and evolution of density perturbations during and 
after chaotic inflation paying attention to the role of the isocurvature mode in a class of 
scalar-tensor theories in which the dilaton is metrically coupled. This paper originates from 
the analysis by Starobinsky and Yokoyama who calculated the spectrum of density 



perturbations from inflation in Brans-Dicke theory. The scalar-tensor theories we study 
include those with negative /3 (the derivative of the Brans-Dicke coupling) which have not 
been studied in the literature ||12[ in which the numerical analysis is limited to positive f3. 
These theories exhibit a significant deviation from general relativity in the strong field region 

i- 

The plan of this paper is as follows: In §2, basic equations of background fields are 
given. Linear perturbation equations are given and solutions to these equations are derived 
in §3. We study the evolution of isocurvature perturbations after inflation in §4. In §5, the 
spectrum of density perturbations are calculated and compared with the COBE observation 
to constrain the parameters of the theory. Finally §6 is devoted to summary. 



2 



II. BACKGROUND FIELD EQUATIONS 



A. scalar-tensor theories of gravitation 

We consider the inflaton, a, coupled to scalar-tensor theories of gravity. The action is 



S 



IGtt 



ah , 



(2.1) 



where ■j/^ is the massless Brans-Dicke dilaton, 'R. is the scalar curvature, is a dimen- 

sionless coupling parameter and V[o^ is the potential of the inflaton. Let us consider a 
conformal transformation, 



9 ah = :;^9ah = e^'^^'^^S'ab, 

with ip = G^^ or a{ilj) = today. Then the action can be written as 



(2.2) 



S 



1 1 1 
irV^-T,9 ^,a'^,h--e 9 cr,acr,6-e V[a) 
Zk^ Z Z 



(2.3) 



where = SnG and if is defined by 



1 fdip 



K J 1p 



We also define 



1 da 1 da 

Kdif Kd(p 

which are related with the PPN parameters {'jEdd, PEdd) and are given as 



(2.4) 



(2.5) 



l7. 



Edd 



l + 2a(v9)2 



Edd 



l + 2a((^)2 



(2.6) 



where the subscript indicates the present value. Here the Brans-Dicke theory corresponds 
to (3{if) = and the general relativity to a{(f) = 0. The recent analyses of the experimental 
data yield |jl4 



I'jEdd - 1| < 2 X 10-^ \4(3Edd - lEdd - 3| < 1 X lO-l 
Combining them yields the following constraint on (3Edd 

|/?£;,,-l|<6xl0-^ 



(2.7) 



(2.8) 
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B. background equations 

We consider the spatially flat FRW spacetime as a background metric 

ds'^ = -dt^ + Ritfdn^. 
The equations of motion derived from the action Eq. (|2.3| ) are 



I) 



K 



0- + 3iJ(T 

H 



(2.9) 



(2.10) 

(2.11) 
(2.12) 

(2.13) 



where ' and ' denote time and a derivatives, respectively. 

The conditions for slow-roll inflation (|o"| ^ |3i^cr|, ^ \3H(f\) are written as 



11-13 



max I -0', -e'V' ) < e^^Via 



1 

2^" '2' 

e-2"V'(a)2 < 6/tV(a)2, 

e-'"\/"((T) < 3KV(a), 

8a{^f < 3, 

|16«(^)2 + 4/?(¥;)| <3. 

When the above inequalities apply, field equations are simplified to 



3H' 



3Hip ~ -4a/te^"V((T) 
3H& ~ -e^'^V'ia). 



(2.14) 

(2.15) 
(2.16) 
(2.17) 
(2.18) 



(2.19) 
(2.20) 
(2.21) 



III. LINEAR PERTURBATIONS 
A. basic equations 

Now we turn to linear perturbations with a perturbed metric in the longitudinal gauge 

ds"^ = -(1 + 2^)dt'^ + R^{1 - 2*)dxl (3.1) 

From the perturbed Einstein equations, each Fourier modeQ satisfies the following equations 
of motion to first order |TT]|T^, 



^ Our convention of the Fourier transformation is ^(k) = (27r) / d^xe ^^^Y{x) for each vari- 
able Y. 
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A;2 



5a + 3H5a + —5a + e^''V"[a)5a = -2e2"y'(a)($ + aK5ip) 
+ (T(4<i> - 2aK5(^ - 2j3K^(p5ip) - 2aKLp5a, 



(3.2) 
(3.3) 
(3.4) 



$ + 4/7$ + (/J + 3if^)$ = y + e^V^fx - e^''V'{a)5a + (e^V^ - 4e^"\/((T))afi:(5¥?J , (3.5) 
together with the Hamiltonian and momentum constraints 



h,2 2 

3H^ +{H + 3H^)^ + — $ = -— \ip5ip + e^^'aSa + e^''V'{a)5a + (e^V^ + Ae^''V{a))aK5^] , (3.6) 

i? 2 L J 



$ + i7$ = y(05v5 + e^V^cr). 



(3.7) 



B. long wavelength perturbations 

Since what we need are the non-decreasing adiabatic and isocurvature modes on large 
scale k <^ RH, which turn out to be weakly time-dependent as will be seen in the final 
result 1^, we may consistently neglect $ and those terms containing two time derivatives. 
Then eqs. (|37^ ), ( p. 41) , and ( p. 71) are simplified to 



1 V 

$ ~ -2aK5ip - --rrSa, (3.8) 

3H6^ ^ -4(3K,'^e^''V6ip, (3.9) 

3H6a - -e^" V6a + 2aKe^^V'5^. (3.10) 
Following [rT|, we find the last two equations can be integrated to give 

^yp^— Qi, (3.11) 

K 

1 v 



5a^--(e-2»Qi + Q2), (3.12) 
hi V 

where use has been made of Eqs.( p.l^ ), ( p.20|) and (|2.21| ). We then find 



^^.Sa'Q^-^(y) (e-^^Qi + g^), (3.13) 



where Qi and Q2 are constants of integration. 
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C. adiabatic and isocurvature modes 



In order to clarify the physical meaning of the above solutions, we should divide them into 
adiabatic and isocurvature modes The primordially isocurvature mode is character- 

ized by its vanishingly small contribution to the curvature perturbation initially, while the 
growing adiabatic mode can be described by the following expressions in the long wavelength 
limit k <^ RH (see Appendix A for a review of its derivation). 









S<Jad 





& 



H 

'r 



R{t')dt' 



H 



(3.14) 
(3.15) 



as a solution of a second-order differential equation Eg. (|A1|) without the source term. Here 
Ci is a constant and the latter approximate equality in each expression is derived by inte- 
gration by parts during the inflationary stage. 

As discussed by Starobinsky and Yokoyama ||rT| when the conditions 

(3.16) 
(3.17) 

hold, one can discriminate between adiabatic and isocurvature modes in the final results by 
defining new constants Ci and C3 as Ci = —Qi — Q2 and C3 = —Q2- Using Eqs.( p.li^ ), 
( ^ ) and we find 


















Ci C3 










5a 




a 


H H ^ 






-2a\ 



(3.18) 

(3.19) 
(3.20) 



In the above expressions, terms in proportion to Ci and C3 represent adiabatic and isocurva- 
ture modes, respectively if the above conditions hold. In more general scalar-tensor theories 
(especially in positive j3 theories) Eq.( |3.16|) and Eq. ( p.lTj ) may not hold and consequently 
terms in proportion to C3 contain the adiabatic mode partially. Then these terms might be 
better referred to as the non-adiabatic mode. But this is a matter of terminology that does 
not affect our final result and we keep using the same expression in this case, too. 



D. quantum fluctuations 

We shall next determine the constants Ci and C3 from amplitudes of quantum fluctu- 
ations of the scalar fields generated during the inflationary stage. Due to the inequalities 
( p.l4| ),( prT5D and ( |2.16| ), Eqs. (|3.3| ) and ( |3.4| ) can be approximated by the equation of motion 
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of a free massless scalar field during inflation for k > RH. The standard quantization gives 
the well-known result, that is, the Fourier components of the fields can be represented in 
the form 



(3.21) 



where tk is the time when /c-mode leaves the Hubble horizon during inflation. Here e^(fc) 
and erj{k) are classical random Gaussian variables with the following averages: 



(e^(k)) = (e.(k)) = 0, (e,(k)e,(k'))=M^^Hk-k'), 2,j = ¥^,a. 
We thus find 



5cj 



-2a jj f^_^' 



- tk 



e-^" , 1 - e 
k + - 



-2a 



-e^(k) 



-3a 



-e.(k) 



-2a 



-e^(k) 



(3.22) 

(3.23) 
(3.24) 



tk 



Note that in Ref. [|T^], they ignored the entropy perturbations at and after the end of inflation 
which corresponds to setting C3 = when they calculated the spectral index although they 
originally took into account these quantities in their derivation. This cannot be generally 
justifled since both 6a and 6ip are stochastic quantities. 



IV. EVOLUTION OF ISOCURVATURE PERTURBATION AFTER INFLATION 

In order to relate the spectrum of $ at the end of inflation to that at decoupling, we 
have to study the evolution of the isocurvature mode together with the adiabatic mode after 
inflation. Unlike the adiabatic modes we do not have a universal formula for the isocurvature 
modes, so we have to calculate them explicitly. 



A. isocurvature perturbation in the reheating stage with Brans-Dicke dilaton 



First let us consider the isocurvature mode during the reheating stage after inflation. In 
this stage, k^V <^ g~2ay// ^g^g Eq.( ^.16D ), so that the inflaton oscillates coherently. Since 

have oscillating coefficients. 



equations of motion of perturbed quantities, ( |3.3| ) through 
we must worry about possible parametric resonance effect. This issue has been studied by 
Kodama and Hamazaki |1^ recently in the case only the inflaton is present and oscillating. 
They have shown that resonance effect is unimportant at least in the long wavelength regime 
which is of our interest. In the present case with two scalar flelds, however, their approach 
is not directly applicable, in particular, for the isocurvature mode |]16[. Hence we have 
numerically solved the evolution equations in the long wavelength limit. Figure 1 illustrates 
an example of the results for the Brans-Dicke theory with uj = 500 and uj = 5. We flnd that 
there is no anomalous growth of perturbation variables and that linear analysis suffices even 
in this regime. We also find qualitatively the same results for more generic scalar-tensor 
theories. 
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B. Isocurvature perturbation in the radiation-dominated stage with Brans-Dicke 

dilaton 



Next let us consider the isocurvature mode during the radiation dominated era after 
the reheating stage. In appendix B, a general expression of the evolution equation for $ is 
given in the dilaton-inflaton-radiation system. After the inflaton has decayed into radiation, 
Eq.( [B37|) becomes 



t.2 

$ + (4 + 3c2)if$ + {2H + 3(1 + cl)H^)<^ + c^— $ 

K 



3 h Vl + Wd 1 + 

where Cg is the sound velocity, 



(4.1) 



P ?>Lp^ + Apr 

with hj = pj + pj, Wj = Pj/pj (j = r,d), and h = p + p. It is now apparent that terms in 
the parenthesis in the right hand side of ( [4. ID are just the entropy perturbation 

y = jA, - -A.. (4.3) 

For the isocurvature perturbation, it is nonvanishing. It is known that if we replace pa with 
the baryon density pb, the isocurvature mode grows as R |T^. Because the dilaton density 
Pd decays as instead of R~^, it is clear that the special solution originating from the 
source term decays as R~^. On the other hand, the homogeneous part of the growing mode 
solution is the same as the adiabatic one. Therefore the isocurvature mode in the radiation 
dominated stage grows no faster than the adiabatic mode. This behavior can be also seen 
by solving the Eq. ( P27|) by using the Green's function method. 

To conclude, the isocurvature mode does not evolve faster than the adiabatic counterpart 
and so is unimportant. The spectrum at the end of inflation is sufficient to compare the 
observation. 



V. SPECTRA OF PERTURBATIONS AND CONSTRAINTS ON PPN 

PARAMETERS 

Now we calculate the density fluctuation spectra produced by inflation. In order to do 
so, we have to specify V{a) and a{ip). We consider chaotic inflation induced by a mass term 
V{(y) = of which we have a sensible particle-physics model in the intermediate scale 

||T8| . Here m is determined from the large-angle microwave anisotropics seen by COBE-DMR 

El. 
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A. general remark on initial conditions 



Initial conditions for Lp and a are given in terms of the number of e-foldings from the 
end of inflation; 

N, = -[ Hdt^^f -^^K^r '^^da, (5.1) 



4 Jip^ a{(p) Ja^ V'{a 

where the suffix e represents the end epoch of inflation. Our present horizon (Hq^ = 
3000/;.~^Mpc) crossed outside the Hubble scale about 60 e-foldings before the end of inflation. 
We assume that the slow-roll condition is satisfied at Nj e-folds from the end of inflation. 



Then the terms containing the two time derivative can be neglected in Eqs. (|3.3| - |3.71 ). Initial 
conditions for H,ip,a are given by Eqs.( |2.19l - p.21| ). 



We need to determine to calculate the spectrum of density perturbations. We nor- 
malize ip in terms of its present value so that ^^(v^) = «o (^^^ ( p.6| )) at present. Then we 
are able to specify v?e by going back to the end of inflation. 

In order to connect with the present length scale, we have to estimate the number of e- 



foldings from the end of inflation to the present pO[. After the end of inflation, the oscillating 
inflaton, which behaves like a dust fluid for the massive inflaton, dominates the universe, 
and then radiation dominated stage follows. The effective temperature, Tj-d, and the scale 
factor /2rd, at the onset of radiation dominated stage is determined from the energy density 

where pe is the energy density at the end of inflation and g^s is the effective number of 
degrees of freedom at temperature T typically taking the value of O(IOO) at Tj-a. Assuming 
the conservation of the entropy in relativistic particles per comoving volume s = ^R^p/T, 
we have 



3 3 3 21 3 3 



43 
IT 



where T^^q is the background neutrino temperature. The number of e-folding from the end 
of inflation to the present is then 

Re 

-72 + iln^-lln^, (5.4) 
3 rripi 3 mpi 

where rripi is the Planck mass. We set the present value of the scale factor Rq to unity for 
convenience. 

Given initial conditions, we solve Eqs.( p.l0| - |2.13| ) numerically using a fourth-order Runge- 
Kutta method. We choose Nf = 65 typically. We express all quantities in units of nipi and 



use X = \nR for the integration variable instead of the cosmic time |20 . 
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The calculations of the perturbed equations Eqs. (|3.3| - |3.7|) are begun when k < RH. The 
initial conditions for S(p and Sa are given by Eq.( 3.2l|) . Since the background fields are in 



their slow-rolling phase, the initial conditions for 5ip, Sa, $, and $ are given by neglecting 
terms with the second order time derivative in Eqs.( p.3|) , (|3^) , (|3.6|) and (p.7|) . We have to 
solve Eq. (|Al|) with or without the source term for the evolution of $ so that we can follow 
the evolution of adiabatic and isocurvature modes separately. The spectral index n defined 
by 

is calculated for the present horizon scale. 



B. COBE normalization 

We have to normalize fluctuations to compare the spectra among different gravitational 
theories. The COBE-DMR data give the normalization. Following the standard treatment, 
the temperature fluctuations are expanded in terms of spherical harmonics as 

AT 

— {e,<p)=J2ai^YiU0,<P), (5.6) 

Im 

and we define 

{al^avm') = ClSll'Smm'- (5.7) 

On COBE-DMR scales, temperature fluctuations are dominated by the Sachs- Wolfe effect 
which is expressed as AT{x)/T = $(x)/3 for scalar perturbations. Here we neglect 
acoustic oscillation terms which may contribute to higher Vs. Therefore Ci is written as 



Ci = (|a/m|^) 

'"^ jfihm), (5^8) 



vri \ 3 ^ 

where r^o = '^Hq^ is the conformal time at present. 

If we approximate the power spectrum of $ by a power law, 

/ ^ S n-l 



fc'(l*r) = ^(;^l , (5.9) 



we can write Eq.(|5.^) for scalar perturbations as 



A r(3-n)r(/+ "-^ 
36 r2(^)r(/ + 

The rms quadrupole Qrms-ps used by the COBE-DMR group is written in terms of C2 as 



Qrms-PS — TqJ (5.11) 
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where Tq = 2.726K. According to the 4 yr COBE-DMR data p 



Qrms-PS = l8fj.K, (5.12) 

for n = 1 Harrison-Zel'dovich spectrum. It takes a shghtly different value if different power 
low indices are considered, but we normalize the spectra by using the above value because the 
index is varied around n = 1 and furthermore the constraints on the scalar-tensor theories 
are imposed primarily from the spectral shape or its power-law index rather than from its 
amplitude. The corresponding amplitude A is 

^ _ 144^ r^(i^)r(/ + ^) f Qr„.s-Ps \' 

^- 5 r(3-n)r(/ + ^) [ To J ■ ^^-'^^ 

Eventually m takes the value ranging from O(10^®mpi) to O(lO^^mpi). 

C. Brans-Dicke theory 

The Brans-Dicke theory corresponds to /3 = and thus 

a(v?) = tto, (5.14) 
a{ip) = aoK(p. (5.15) 

Initial conditions for and a are 

n{p = iaoN + nipe, (5.16) 
K(pe ~ lOao, (5-17) 



3 20^0 



In Figs. 2, the spectra for each ao or uj are shown. In the Brans-Dicke theory, 

the conditions Eqs. (|2.14|) and (|2.15|) hold and therefore we are justified to say that terms in 



proportion to C3 represent the isocurvature mode. We see that considering the current limit 
on uj {uj > 500) in Eq.( p.l8 ), the contribution of the isocurvature mode is totally negligible 



in the Brans-Dicke theory in agreement with Starobinsky and Yokoyama . 

The spectral indices on the comoving horizon scale today are n = 0.962, 0.966, 0.967 for 
u = 500, 5000, 50000, respectively. Note that n = 0.967 in general relativity. 



D. a class of scalar-tensor theory 

We take the following simple functional form for a{(f) in accord with Damour and 
Nordtvedt [0] 



= ^i^'^iv'^ - V^l), ^ = const, (5.19) 
a{ip) = PKif. (5.20) 
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We choose the origin of the field ip so that a{(p) = at present in order to reproduce the 
correct value of the gravitational constant .0 We have to determine the value of p at the end 
of infiation, i.e., ipe-, to calculate the spectrum of density perturbation. Since we normalize 
ip in terms of its present value, we need to go back to the end of infiation to specify <^e- The 
evolution of ip in the radiation or matter-dominated universe is analyzed by Damour and 
|T^(see also ||2^). In the radiation-dominated stage (p is hardly changed. 



Nordtvedt |T^(see also 12^). In the radiation-dominated stage (p is hardly changed. We 
assume for simplicity that ip^ is equal to ip at the beginning of matter dominated stage. For 
<C l(see Eq.( p.l8| )), (pe is determined by 



H(pe 



/3 ■ 



Then (p is given in terms of as 



Kip ~ Kipf.e ' ~ — e 



4/3Ar+10/3 



(5.21) 



(5.22) 



For a with (see Eq.(g7[5|)) 




(5.23) 



o is given by 



3 2(3 



(5.24) 



where Ei{x) is the exponential integral function defined by Ei{—x) = — e~^dt/t. 

In Figs. 3, the spectra for several /3 with ao = 10~^ are shown. In some cases the mode 
coming from C3 part in Eq.( |3.18| ) is not negligible and the spectrum can be different from 
fiat one. In scalar-tensor theories with positive (3, the conditions Eqs. (|3.16| ) and (|3.17|) do 
not hold since initially a can be greatly different from zero (see Eq. (|5.22| )), and terms in 
proportion to C3 contain the adiabatic mode in addition to the isocurvature mode. For 
theories with negative (3, however, terms in proportion to C3 represent the isocurvature 
mode. 

In Fig. 4 we show that contour plot of n on ~ plane. We see that for /? > n 



falls down very steeply with increasing /3. Thus the observed spectrum by COBE-DMR [|T9 

n 



1.2 ± 0.3, strongly constrains (3. This possibility was pointed out in |T2| for /3 > 0. 
Interestingly, the contour level can be fitted by a linear function in logao- For n > 0.9, 



9.0 X 10"^ log 



"0 



V5 X 10-2 



- 1.5 X 10"^ < P < -8.9 X 10"^ log 



ao 



5.3 X 10~l 
(5.25) 



2 The authors of [12| take a{(p) = oi + a2K(f> and choose the origin so that a = oi at the end of 
inflation. However, a deviates from oi afterwards and may conflict with the limit on oq. 
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For n > 0.7, 



5.0 X IQ-^ log ( I - 3.7 X 10^2 <f3< -9.2 x 10"^ log ( ) + 7.3 x 10^1 

Vv5xl0~2y y V5 X 10-2/ 

(5.26) 

Since the tensor mode of perturbations does not contribute to the spectrum significantly 
in the above range of the spectral index, the above inequality gives the constraint on f5. 
Note that the constraint on (3 has little dependence on ao as well as on the details of 
normalization of fiuctuations; adopting the inflaton mass which is ten times larger results in 
only a few percent changes in the constraint. We also note that the constraint is stronger by 
an order-of- magnitude than those obtained by post-Newtonian experiments(see Eq.( p.8| )). 



VI. SUMMARY 

We have studied density perturbations produced during chaotic infiation taking both 
adiabatic and isocurvature modes into account in a class of scalar-tensor theories of gravity 
in which the dilaton coupling is metric one. The spectrum of the density perturbation 
produced by chaotic infiation in a scalar-tensor theory can be a non-fiat one because of the 
variable Brans-Dicke coupling. The spectrum observed by COBE-DMR thus constrains such 
a coupling. Assuming a simple coupling function, the constraint is found to be stronger by an 
order-of-magnitude than those obtained by post-Newtonian experiments. It is interesting 
to note that in the case of the simple coupling function a{ip) = jSip employed here, our 
constraint is much more stringent than that obtained by the binary pulsar experiment ||23|| , 
/3 > —5, in a sense. Indeed ao must be ao < 10~^°°°(!) so that the constraint Eq. (|5.25|) or 



Eq.( [5.26D is looser than (3 > —5 if these fitting formulas are extrapolated thus far. Thus we 



may conclude that cosmological infiation is suited in the framework of general relativity, or 
more advocatively, infiation favors Einstein. 
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APPENDIX A: FORMULAE FOR ADIABATIC MODE 



In this appendix we review derivation of the formulae for the adiabatic mode Eqs.( p.l4D 
and ( 3.151 ). See @] for a different approach. 



Using Eqs.(pl), (^) and (^), we get 
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$ + (4 + ?,cl)H^ + {2H + 3(1 + cl)H^)^ + c^— $ 

.2 1 N ^ , ,„2 ^ nrl- A -la. 



- 1)^'^' + ^ -2^ 2a-2 (^V - 4e^'^««:\/a)(a5v^ - 05^), (Al) 



i?2 (^2 + g2a^ 

where is defined by 

2e^"(4«/.\/v9 + ra) 
^"p" + 3ff(02 + e2V2) • ^^^^ 

The left-hand-side of Eq. (|A1|) is just the same as that of the adiabatic perturbation for hydro- 
dynamical matter and in case of the adiabatic perturbation right- hand- side term vanishes. 
For long wavelength perturbations right-hand-side vanishes when 

^ = ^. (A3) 

a 

To derive the solution for the adiabatic mode Eq.( |3.14] ), it is convenient to use the constancy 
of the Bardeen's C 

C = ^($ + H-^^) + $ = const = Ci, (A4) 

H 

which follows from the fact that the left-hand-side of Eq. (|AT|) is equal to —H(/H in the 
long-wavelength limit. Substituting $ — Ci = the above equation becomes a homogeneous 
form 

$ + iJ<l - — $ = -CiH. (A5) 
H 



"(^i), (A6) 



Noting the left-hand-side is written as 

'r \H 

the solution is immediately given as 

$ = -Ci^ J R{t')dt' + C2^, (A7) 

where C2 corresponds to the amplitude of the decaying mode which is neglected in Eq. ( P.14| ). 
To derive Eq.( |3.15D , we use Eq.(p.7|) as 



$ + if$ = -Ci— / Rdt' + C2 
R J 



{^5^ + e^VH = -^(02 + e^V^). (A8) 



Employing Eq.( p.l3| ) we find 



Lf) a R 
Neglecting the decaying mode, we get Eq. (|3.15|) . 



^-^ = - = \,[C,lRdt'-C2]. (A9) 
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APPENDIX B: LINEAR PERTURBATIONS FOR 
DILATON-INFLATON-RADIATION SYSTEM 



Basic equations in the Einstein frame are the Einstein equation 
Gab = «^ (Tab + da(fidb(fi - ^^QabQ^'^dc^dd^ + e^"9„(796(7 - gabi^g'^e^'^dcadaa 

the dilaton equation of motion 

Dip = -Ka{^){T - g^'^e^^'dcadaa + 4e^"l^(a)), 
the matter equation of motion (Bianchi identity) 

W,T\ = Ka{T - g^'^e^'^dcadda + 4e^«F(a))Va<^, 
and the inflaton equation of motion 

where — SttG. 

1. background equations 

The matter energy momentum tensor is given by 

Tab = PrUaUb + Pr{9ab + UaUb). 

Assuming the flat FRW model, the background equations are 

,2 



2 

(fi + 3Hip = Ka{-pr + 3pr + e^"(7^ - 4e^"y((7)), 

{PrR^) + Pr{R')' = {pr - Zpr) a{ip) , 

where Pr = Pr/'i- 

2. linear perturbations 

We employ the metric perturbation in the longitudinal gauge as 
ds^ = -(1 + 2^)de + R^{1 - 2$)(ix2. 
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Then the perturbed Einstein tensors are [p4 



3R^^ + 3Rm + 



m + m 



{R^ + 2RR)^ + Rm + {R^^y + Rm 

As for the matter components, hydrodynamical perturbations of radiation are 

6T\ = -pA 



6T 



{Pr+Pr)vYj, 



The dilaton perturbations are expressed as 

k 



6r, 

and the inflaton perturbations as 



6T\ 



6T' 



5T 



k 



R 



Thus we have the perturbed Einstein equations 



5!- 



R 

(pr + Pr)v— + (pSip + e^"'&6a 
k 



The perturbed equations are rewritten as 
$ + (4 + 3cl)H^ + (2i7 + 3(1 + cl)H^)^ + c^— $ 



+ 2e^"\/'((j)5(T + 8/€ae*'^\/(a)5v3 



(B12) 
(B13) 
(B14) 

(B15) 
(B16) 
(B17) 

(B18) 
(B19) 
(B20) 

(B21) 
(B22) 

(B23) 



(B24) 
(B25) 

(B26) 



, (B27) 
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where is defined by 



(B28) 



We have arranged the left-hand-side of Eg. ( |B27] ) so that it is just the same as that of the 
adiabatic perturbation for hydrodynamical matter. 

Now let us examine the structure of the right-hand-side of Eq. (|B27|) in detail. Before 
doing so, we define the following useful gauge invariant variables [Q. For a general multi- 
component system {pa,Pa), and Va are defined by 



= -pa5a, 

6T^j = haVaYj, 



(B29) 
(B30) 



where ha = Pa + Va- The associated gauge invariant variables Aq and Va are defined by 



/\a = 5a + 3{l+Wa){l-qa) 
Va = Va, 



-Vt 



(B31) 
(B32) 



where Wa = Pa/ Pa and Vt = T,haVa/h with h = T,ha and qa is the source term for each 
energy momentum tensor. For the case we are considering, i.e., the radiation- inflaton-dilaton 
system, the respective Aa(Ar, Aj, A^) are written as 



A, 



1 + Wr 

A, 

1 + Wi 

Arf 



4 k 



1 + 



4e^"\/(a))Ka0' 



3Hhr 



Vt, 



6a 
& 



RH, 



$ H h Ka5ip + 3-PVt, 

k 



:$ + — + 3^ 

l + Wd ^ k 

Then Eg. ( [B27| ) can be rewritten as 



3Hh, 



Vt 



(B33) 
(B34) 
(B35) 



^2 

$ + (4 + 3cl)H^ + {2H + 3(1 + c^,)H^)^ + c^— $ 

R'^ 



id - r 



R^ 



2 r 1 /? 

+ 2hrVT 



+ 2^e2"\/'(a)a 



6a vR^ 
a k 



+ 



h 



Kae^''V{a)ip 



'Sip 



h \aip 



vR" 



(B36) 



When radiation is absent the expression agrees with that of inflaton-dilaton system. 
By expanding out the term (c^ — we finally have 
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t.2 

$ + (4 + 3c^)i/$ + i2H + 3(1 + c^)i/2)$ + $ 



2 l3/i 



A, 



1 + w,. 1 + 



+ hrhd 



A, 



A„ 



1 + tVr 1 + 



2RVt^^2 



3(7 A; / /i \3<^ A; 



h 



(B37) 
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FIGURE CAPTION 



Fig.l. Evolution of the inflaton and the curvature perturbations during the oscillation regime 
of the inflaton in the Brans-Dicke theory with {a)uj = 500 and (b)a; = 5. The abscissa 
is the e-folding number after the end of inflation. 

Fig. 2. The spectra of perturbations in Brans-Dicke theory with u = 500, 5000, 50000. The 
abscissa is the wavenumber in unit of /iMpc^^. Logarithm is base 10. The isocurvature 
mode is negligible in this theory. 

Fig. 3. The spectra of perturbations in scalar-tensor theories with = 1-0 x 10^^ and 
P = 0.02,0.005,-0.04. For larger the mode coming from C3 part in Eg. ( p. IS]) 
is not negligible and the spectrum has slope. The respective spectral index is n = 
0.640,0.967,0.734. 

Fig. 4. The contour plot of spectral index n for scalar-tensor theories. The contour levels 
n = 0.9, 0.7 are shown. The region bounded by these courves is the allowed region. 
We find non-zero f3 theories are strongly constrained. The shaded region is excluded 
by the solar system experiments; the vertical line is the constraint by the expetiment 
of the deflection of light (or the time delay of light): ao < V5 x 10~^. The dotted 
curve in the shaded region is the constraint by the lunar-raser-ranging experiment: 
al\/3o\ < 3 X 10"^ 
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